Optical Conductivity a(u) and Resistivity pd c 
of a Hole Doped Spin-Fermion Model for Cuprates 



Mohammad Moraghebi 1 , Seiji Yunoki 2 and Adriana Moreo 1 
1 Department of Physics, National High Magnetic Field Lab and MARTECH, 
Florida State University, Tallahassee, FL 32306, USA 
2 Istituto Nazionale di Fisica della Materia and International School for Advanced Studies (SISSA), 

via Beirut J^, Trieste, Italy 
(February 1, 2008) 

The optical conductivity and Drude weight of a Spin-Fermion model for cuprates are studied as 
a function of electronic density and temperature. This model develops stripes and robust D-wave 
pairing correlations upon hole doping, and it has the advantage that it can be numerically simulated 
without sign problems. Both static and dynamical information can be obtained. In this work, it 
was possible to analyze up to 12 x 12 site clusters at low temperatures ranging between O.Olt and 
OAt (between 50K and 500K for a hopping amplitude t~0.5 eV). As the temperature is reduced, 
spectral weight is transferred from high to low frequencies in agreement with the behavior observed 
experimentally. Varying the hole density, the Drude weight has a maximum at the optimal doping for 
the model, i.e., at the density where the pairing correlations are stronger. It was also observed that 
the inverse of the Drude weight, roughly proportional to the resistivity, decreases linearly with the 
temperature at optimal doping, and it is abruptly reduced when robust pairing correlations develop 
upon further reducing the temperature. The behavior and general form of the optical conductivity 
are found to be in good agreement with experimental results for the cuprates. Our results also 
establish the Spin-Fermion model for cuprates as a reasonable alternative to the t — J model, which 
is much more difficult to study accurately. 

PACS numbers: 74.25.Gz, 71.10.Fd 



I. INTRODUCTION 

Infrared measurements are an important probe of the 
dynamical properties of the high critical temperature su- 
perconducting cuprates. In particular, the real part of 
the optical conductivity <j(uj) provides useful insight into 
the electronic structure of these materials. Measure- 
ments of the optical conductivity have been performed in 
both hole and electron doped cuprates. Common features 
observed include the absence of absorption below the 
charge transfer gap in the insulating phase (half- filling). 
It has also been reported the rapid transfer of spectral 
weight to low frequencies with increasing hole doping, 
which gives rise to a Drude-like peak at lu — and 
a broad mid-infrared feature, while the spectral weight 
above the charge transfer gap decreases. |||| The re- 
ported integrated conductivity up to 4 eV (beyond the 
charge transfer gap) in La2- x Sr x Cu04 remains approx- 
imately constant with doping, indicating that spectral 
weight is redistributed from the charge transfer band to 
lower frequencies. H It is also observed that at low fre- 
quencies <j{oj) ~ instead of the 1/ui 2 behavior ex- 
pected in a standard metal. 

A number of the above mentioned features are well 
reproduced by models of strongly correlated electrons, 
including the Hubbard, t — J, and related Hamiltoni- 
ans, (JJ^] even in regimes where superconductivity was 
not numerically detected in the ground state. This sug- 



gests that some of the previously described properties 
of cuprates may just be the effect of strong correla- 
tions among the electrons, rather than superconductiv- 
ity which is expected to develop at much lower tem- 
peratures that those previously studied. A problem in 
the calculation of the optical conductivity in Hubbard 
and t — J models is that only small clusters of about 
20 sites can be studied with exact diagonalization tech- 
niques. KM Alternative approaches, such as the quantum 
Monte Carlo method allows the analysis of larger clus- 
ters (8x8), but at very high temperatures and with less 
precision. However, recent investigations have shown 
that these problems can be overcome by using a Spin- 
Fermion model (SFM), which is much easier to study 
numerically. In fact those previous investigations have 
reported the presence of a stable striped ground-state 
upon doping, and robust D-wave pairing correlations in 
the SFM. [f|-[jTj The presence of phenomenologically ob- 
served regimes of the high-T c cuprates in a simple model 
of interacting mobile carriers and localized spins is re- 
markable, and opens the way to detailed numerical inves- 
tigations of relevant physical quantities. In this context 
clusters as large as 12 x 12 can be investigated for low 
temperatures ranging between O.Olt and O.lt, consider- 
ably improving upon size and temperature limitations of 
the t — J and Hubbard models. In particular, the ef- 
fect of the stripes on a(uj) and its distribution of spectral 
weight, and dependence with temperature can be studied 
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in this case, as reported here. 

The organization of the paper is as follows: In sec- 
tion II the Hamiltonian and the numerical technique are 
described. Results at T « are discussed in section 
III, while section IV is devoted to the finite temperature 
analysis. An estimation of the linear dc resistivity is pre- 
sented in section V and the conclusions appear in section 
VI. 



II. MODEL AND TECHNIQUE 

The SFM is constructed as an interacting system of 
electrons and spins, which mimics phenomenologically 
the coexistence of charge and spin degrees of freedom 
in the cuprates. (IsJlJB. Its Hamiltonian is given by 
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where cj creates an electron at site i = (i x ,i y ) with 



< -^.ff-ypCip is the spin of the 



spin projection 7, s ; =Y] 7 ^ ^ 
mobile electron, the Pauli matrices are denoted by er, 
Si is the localized spin at site i, (ij) denotes nearest- 
neighbor (NN) lattice sites, t is the NN-hopping ampli- 
tude for the electrons, J > is an AF coupling between 
the spins of the mobile and localized degrees of freedom, 
and J' > is a direct AF coupling between the local- 
ized spins. The density (n) = l — x of itinerant electrons 
is controlled by a chemical potential ji. Hereafter t = 1 
will be used as the unit of energy. From previous phe- 
nomenological analysis, and as in previous papers p,^2|, 
the coupling J=2 and the Heisenberg coupling J'=0.05 
are selected, since they provide the striped and pairing 
states observed in cuprates. To simplify the numerical 
calculations, avoiding the sign problem, localized spins 
are assumed to be classical (with |5i|=l). This approx- 
imation is not as drastic as it appears, and it has been 
extensively discussed in detail in previous publications 
such as Ref. The model will be studied using a 

Monte Carlo method, details of which can be found in 
Ref. p5|. Periodic boundary conditions (PBC) are used. 
The real part of the optical conductivity is calculated 

as ira 
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where a = x, y; Z is the partition function; |n) and |m) 
denote eigenstates of the system; and J a is the current 
operator given by 



(3) 



Since we are using PBC, the Drude weight D a is calcu- 
lated indirectly, using the two-dimensional sum rule, 
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where (—K a ) is the average kinetic energy in the direc- 
tion a and V is the number of sites in the lattice. Assum- 
ing the existence of a contribution DS(uj) at zero energy, 
we obtain 



7T 



(-K a ) 1 



2N 



0+ 



du>cr aa {uj), 



(5) 



This procedure to calculate the Drude weight is common 
ractice and it has been discussed extensively in reviews. 
In this work, both directions x and y must be analyzed 
independently since at low temperature, T < 0.025i, the 
symmetry under rotations is broken due to the formation 
of stripes. As a consequence, the optical conductivity will 
be measured along the direction parallel and perpendic- 
ular to those stripes. 

In numerical calculations of the optical conductivity 
it is customary to replace the delta functions in Eq.(2) 
by Lorentzians of width e. The smearing of the deltas 
helps to simulate effects not considered in the microscopic 
Hamiltonians, such as disorder which contributes to the 
experimentally observed broadening of the peaks. This 
peak broadening also occurs naturally in the Monte Carlo 
simulations since the position of the poles varies slightly 
between successive iterations. Here, the incoherent part 
of the optical conductivity is calculated by effectively giv- 
ing a width e w O.Oli to the delta functions. That same 
value of e will be used in the Lorentzian for the Drude 
weight. 



III. RESULTS AT T~0 
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FIG. 1. (a) The optical conductivity versus uo at T = 0.01 
in the direction perpendicular to the stripes and different den- 
sities on a 12 x 12 cluster, (b) Same as (a) but in the direction 
parallel to the stripes. 
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In Fig. la, we present the real part of the optical con- 
ductivity measured in the direction perpendicular to the 
stripes on a 12 x 12 lattice at T — O.Oli for different values 
of the electronic density (n) . At half-filling all the spec- 
tral weight appears beyond lu w 3.8i. This corresponds 
to about 2 eV if we assume t w 0.5eV, in agreement with 
experimental results. This half- filling gap is created by 
the coupling J, which plays a role analogous to U in the 
Hubbard model. In fact, J suppresses double occupancy 
of the mobile carriers as strongly as U does. This impor- 
tant point has been extensively discussed particularly in 
the manganite literature where a similar model, but with 
different signs and values of couplings is used to describe 
e g electrons moving in a t 2g background. [jl7| 




<n> <n> 



FIG. 2. (a) The Drude weight as a function of the electronic 
density at T = 0.01 for currents perpendicular (parallel) to 
the stripes denoted by filled (open) circles, (b) The kinetic 
energy as a function of the electronic density at T = 0.01 in 
the direction perpendicular to the stripes (filled circles) and 
parallel to them (open circles) . (c) Integral of the incoherent 
spectral weight of <j(lu) versus the electronic density, calcu- 
lated along the direction perpendicular to the stripes (filled 
circles) and parallel (open circles). 

For (n) = 0.875, spectral weight is observed at low 
values of to, again as in the experimental data for the 
cuprates. In the same figure, data for (n) = 0.80, 0.75 
and 0.68 is shown. It can be observed that spectral 
weight is transferred from above the gap to low ener- 
gies in all cases. At the same time, the Drude weight, D, 
shown in Fig. 2a with filled circles, becomes finite with 
doping, indicating that the system transforms from an 
insulator to a conductor. Notice that D reaches its max- 
imum value for (n) w 0.75 — 0.80. According to previ- 
ous investigations, this corresponds to the optimal dop- 
ing for the SFM, i.e., these are the densities at which the 
D-wave pairing correlations are the strongest. |ll| The 
maximum obtained in the Drude weight is due to the 
interplay between the kinetic energy and the incoherent 
spectral weight, which are shown in Fig. 2b and c respec- 
tively. It is interesting to notice that in previous studies 
performed on the large U Hubbard and t — J models, it 
was observed that the kinetic energy reaches a maximum 



at quarter-filling in regimes without stripes. |18(] How- 
ever, in the SFM, the kinetic energy perpendicular to the 
stripes has a maximum for (n) sa 0.7 (see Fig. 2b). This is 
in agreement with experimental results which found that 
the plasma frequency, proportional to the kinetic energy, 
grows with doping in the underdoped region and stops 
growing at optimal doping. |lj| This also establishes an 
important difference between previous t—J model studies 
and those reported here. 

It has also been observed that in the direction parallel 
to the stripes the conduction is much suppressed, as it 
can be seen in Fig. lb. This is in agreement with previous 
dynamical studies of the model, which showed that the 
kinetic energy is larger in the direction perpendicular to 
the stripes flQ] , as appears also in Fig. 2b. It is also clear 
that the kinetic energy continuously decreases with dop- 
ing in the parallel direction. The corresponding Drude 
weight is denoted with open circles in Fig. 2a. 

IV. RESULTS AT FINITE TEMPERATURE 

The next step is to study the dependence of the optical 
conductivity with temperature. Experimental measure- 
ments have been carried out mainly in a narrow range 
between 0.01 and 0.25 eV, which corresponds tow < 0.50 
in our scale. For La2_ x Sr x Cu04 (LSCO) measurements 
have been reported for x=0.13 and 0.14 (underdoped) 
and x=0.22 (overdoped) at temperatures ranging be- 
tween IQK and 400K, @ while underdoped Bi2212 and 
optimally doped Y123 have also been studied. pl[ 
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FIG. 3. (a)The optical conductivity of the Spin-Fermion 
model versus lu at (n) = 0.75 in the direction perpendicular 
to the stripes and different temperatures on a 12 x 12 cluster. 
The temperatures are indicated; (b) Detail of the curve shown 
in part (a) for T = O.Olt where the data for 0.1 < lu < 0.25 
are fitted by an 1/lu curve. 

In the underdoped case, it has been experimentally 
observed that as the temperature decreases from A00K, 
spectral weight is transferred from intermediate toward 
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lower frequencies. This depletion of spectral weight is 
associated with the opening of a pseudogap in the den- 
sity of states at a temperature T* estimated to be above 
400AT. In the overdoped and optimally doped samples, 
on the other hand, the depletion of spectral weight is 
observed below T c indicating that T* « T c in this case. 

In Fig. 3a, the optical conductivity versus u> is shown 
at optimal doping (n}=0.75 for different temperatures. 
Based on the previous studies of D-wave pairing cor- 
relations, the critical temperature is T c «0.025. It can 
be clearly seen that spectral weight is transferred to 
lower frequencies as the temperature decreases. The mid- 
infrared (MIR) weight increases, as well as the Drude 
weight, whose inverse is displayed in Fig. 4 as a function 
of temperature. 

The MIR feature is characteristic of doped cuprates 
and it is located atw« 0.3cV. This is in very good 
agreement with our results presented in Fig. la and 3a 
where the MIR feature appears for oj w 0.5 — 0.8, which 
corresponds to 0.25-0.40 eV for t — 0.5eV. 

As mentioned before, it was experimentally observed 
that at low frequencies (above ui « 0.03eV) the decrease 
in the conductivity is closer to l/u> than 1/uj 2 , which 
would have been expected for free carriers. The best fit 
of our data for ct(cj) in the range 0.1 < u> < 0.25 with a 
power law C/w Q corresponds to 0.9 < a < 1.3, in all the 
range of temperatures and dopings studied, in excellent 
agreement with the experimental results. In Fig. 3b we 
show, as an example, the 1/lu fit for the data at optimal 
doping and T = 0.01£. It appears that the mixture of 
Drude weight and MIR band at low frequencies conspires 
to produce the 1/uj behavior. 



V. ESTIMATION OF THE RESISTIVITY 



termine the onset of superconductivity by studying l/D 
as a function of the temperature since we cannot dis- 
tinguish between a perfect metal and a superconductor 
by monitoring the Drude weight in this case. However, 
we observe a rapid reduction of l/D at T c (see Fig. 4). 
This reduction occurs as the D-wave pairing correlations 
develop a robust value at long distances. The longest 
distances in our clusters is L=6 and as a consequence 
D(6) is also shown in Fig. 4 as a function of the temper- 
ature (open circles). In this case T c = 0.025£ « 150AT, 
which is the temperature for which a pseudogap opens in 
the density of states, completing an impresive agreement 
with cuprate's phenomenology. 
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FIG. 4. The inverse of the Drude weight (filled circles, scale 
on the left) as a function of the temperature, for the opti- 
mal doping (n) — 0.75, and for currents perpendicular to the 
stripes. The line indicates a linear fit. The open circles indi- 
cate the D-wave pairing correlation at the maximum distance 
D(6) (scale on the right). The qualitative agreement with 
cuprate's experiments is evident. 



It has already been mentioned that the effects of dis- 
sipative processes that would create a finite resistance in 
a metal have to be added ad-hoc, by replacing the delta- 
functions at lu — by Lorentzians, since these processes 
are not included in the Hamiltonian of the Spin-Fermion 
model. Once this is done, the inverse of the Drude weight 
should provide information about the dc resistivity pdc- 
|^| Transport experiments indicate that a linear behav- 
ior occurs for T > T* , while in the pseudogap region 
of the underdoped compounds a faster decrease of pdc is 
found for T c < T < T* . @|1 A power law dependence 
Pdc ~ T 1+s was observed in the overdoped regime with 
S = 0.5 at 34% doping in LSCO which is optimally doped 
at 15%. H 

It is remarkable that many of the above mentioned ex- 
perimental characteristics are qualitatively captured by 
the SFM. In particular, in Fig. 4, a linear fit of our data 
at optimal doping is obtained above T c (filled circles). 
Notice that, in principle, since no dissipative processes 
are included in the Hamiltonian it is not possible to de- 



Behavior in qualitative agreement with the experi- 
ments is also observed away from optimal doping. In 
our underdoped regime, with a representative density 
(n) = 0.875, the inverse of the Drude weight and the 
D(6) pairing correlation are shown as a function of tem- 
perature in Fig. 5a. In this case, linear behavior reducing 
T is observed up to T = 0.033i «s 200A" when a pseu- 
dogap starts developing in the density of states. Thus, 
this temperature may be associated to the experimental 
T*. Below this temperature, the slope of l/D increases 
and a further reduction of the resisitivity is observed at 
T c = 0.017i « 100A which is the temperature where 
the D(6) pairing correlations start to develop. This is in 
qualitative, although not quantitative, agreement with 
the experimental data for La2- x Sr x Cu04. 

Data in the overdoped regime, for (n) — 0.68, are dis- 
played in Fig. 5b. As in the previous two cases, a sharp 
decrease in l/D is observed at T c when the D(6) pair- 
ing correlations start to increase. However, above T c , we 
observe indications of superlinear behavior like experi- 



4 



mentally observed 23 1 , particularly for the points at the 
lower temperatures. 
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FIG. 5. (a) Same as Fig.4 but for (n) = 0.875 (under- 
doped). The continuous line is a linear fit of the 1/D points 
above T*; (b) Same as Fig.4 but for (n) = 0.68 (overdoped). 
The continuous line is a T 1 ' 5 fit of the four first 1/D points 
above T c . 

On the other hand, while the onset of superconductiv- 
ity is very obvious in resistivity measurements, it is not 
clear that the incoherent part of the optical conductivity 
is sensitive to it, particularly in the overdoped regime. 
The transference of spectral weight from higher to lower 
frequencies observed as the temperature decreases has 
been associated to the opening of the pseudogap rather 
than to the onset of superconductivity. No particular 
feature is observed at T c in this case. [g0| 

An interesting characteristic of the SFM results is that 
spectral weight is transferred to lower frequencies when 
the temperature decreases, as experimentally observed in 
the cuprates. This behavior is not obvious. In fact, since 
the SFM develops stripes, i.e., charge ordering, at low 
temperatures one possible behavior would have been that 
spectral weight was transferred toward higher frequencies 
when the temperature decreases which is the expected 
behavior for systems that develop charge density waves. 



VI. CONCLUSIONS 

Summarizing, the optical conductivity of the SFM ap- 
pears to have many features in common with those of 
the high T c cuprates. A transference of spectral weight 
from high to low frequencies is observed upon doping. 
The Drude weight reaches its maximum value at optimal 
doping and its inverse, roughly proportional to the resis- 
tivity, decreases linearly with temperature in this regime. 
Qualitative agreement with the experimental behavior of 
the resistivity is also observed in the underdoped and 
overdoped cases. The simplicity of the SFM allows the 
study of larger clusters than with the Hubbard and t — J 
Hamiltonians, and a much broader range of tempera- 
tures can be explored as well. Our results establish the 



Spin-Fermion model as a qualitatively realistic model for 
cuprates. 
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